Instabilities in droplets spreading on gels 
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We report a novel surface-tension driven instability observed for droplets spreading on a compliant 
substrate. When a droplet is released on the surface of an agar gel, it forms arms/cracks when the 
ratio of surface tension gradient to gel strength is sufficiently large. We explore a range of gel 
strengths and droplet surface tensions and find that the onset of the instability and the number of 
arms depend on the ratio of surface tension to gel strength. However, the arm length grows with 
an apparently universal law L oc t 3 ^ 4 . 
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The surface-tension driven spreading of liquids is in- 
dustrially and biologically important, and has been stud- 
ied in detail on both solid and liquid substrates [J, H, 0| • 
Less is known about how droplet spreading is modified in 
the presence of a compliant substrate, a situation espe- 
cially relevant to biological applications [4, 5]. We per- 
form droplet-spreading experiments on gel agar, a vis- 
coelastic material, to explore the influence of substrate 
on the spreading dynamics of the droplet. We find a 
novel branching instability with an onset that is con- 
trolled by the ratio of surface tension difference to the 
shear strength of the gel. The existence of a spreading 
morphology in which a spreading droplet becomes spa- 
tially localized has important implications for the indus- 
trial and medical application of surfactants. 





(b) ; 

J 

^ 5mm 




\ 5 mm 1 






Droplets spread differently on liquids, which are mo- 
bile, than on solids, which are essentially rigid [1, @, 0|- 
The present experiments on spreading on a viscoelastic 
substrate (gel agar) are intended as a way to span these 
two limits: by increasing the agar concentration, one can 
tune the substrate from liquid-like to solid-like behav- 
ior. The difference in surface tension between the droplet 
(PDMS, a silicone oil; or Triton X-305, a surfactant so- 
lution) and the substrate drives the spreading process. 

Previous work on the spreading of droplets on gel sub- 
strates [1, Q showed circular spreading with rates in- 
termediate between those observed on solid and liquid 
substrates, contrary to what we find. In addition, prior 
studies of viscoelastic substrates have focused on sub- 
strate breakup or fracture when subjected to stresses 
firl Ull . [lil Ha |. Here, we observe a failure of the gel 
driven by surface tension. After the initial failure, there 
are two morphologically different manifestations of the 
instability, influenced by both the substrate elasticity 
and the surface tensions, which we call starbursts and 
wispy drops. For weak gels (shear modulus G < 30 Pa), 
the drop breaks into distinct crack-like spreading arms 
in a starburst formation, as shown in Fig. [T^i. Morpho- 
logically, this is similar to cracking patterns observed in 
[loL [l2| . We observe that the arms have steep sides ex- 
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FIG. 1: Sample images of droplets spreading on gel agar sub- 
strates, (a-c) Shadowgraph images adapted from [14|: bright 
lines represent the outer boundaries of the arms, (a) Triton 
X-305 solution starburst: droplet concentration cj> = 50 ppm 
(a = 63 mN/m), X = 0.12%w (G = 16 Pa), and V = 5/jL 
after 14 sec. (b) Triton X-305 solution wispy drop: droplet 
concentration 4> = 100 ppm (a = 60 mN/m), X = 0.14%w 
(G = 28 Pa), and V = 5fiL after 27 sec. (c) PDMS starburst: 
a = 20 mN/m, X = 0.08%w (G = 1.3 Pa), v = 1000 cS and 
V — 5pL after 5 sec. (d) Scattering image of Triton X-305 
starburst (inset) and cross section of arm. 



tending into the gel and are filled with material from the 
spreading droplet. Above the onset of the starburst insta- 
bility, the rate of spreading is found to be controlled only 
by the width of the arms, with collapse in the data across 
substrate modulus, surface tension, viscosity, and droplet 
volume. Well above the gelation transition (G > 30 Pa), 
the drops are typically circular as on a pre- wet solid sub- 
strate, i.e. there is no indication of failure. However, 
after long times some droplets (on gels near G » 30 
Pa) demonstrate the wispy morphology shown in Fig.[T|3. 
We investigate the morphology and dynamics of the star- 
bursts in terms of droplet surface tension, viscosity, and 
volume, and the strength of the gel. 

The substrate is made from agar, a polymer composed 
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of subunits of galactose (a sugar); agar undergoes a phys- 
ical gel transition for weight concentrations of agar above 
0.014% weight in water at 20.0°C [I3|. We express agar 
concentrations as X = M agar /(M agar + M water ) %w, 
where M is the mass of each component of the gel. We 
examine droplets spreading on rs 5 mm thick gels with 
concentrations 0.06 < X < 0.14. These values corre- 
spond to a shear modulus of 1 Pa < G < 30 Pa [13] and 
a surface tension 67 mN/m < a < 75 mN/m. We mea- 
sure the latter using the pendant drop method [3]. For 
X < 0.06%w the droplets spread out as on water, and 
for X > 0.014%w they remain circular and small (mm 
radius) [13]. 

We use two types of droplets, Triton X-305 (octylphe- 
noxy polyethoxy ethanol from Dow Chemical) at con- 
centrations 4> = 1 to 1000 ppm in deionized water and 
pure polydimethylsiloxane (PDMS) at several viscosities 
and surface tensions. Triton X-305 is a non-ionic surfac- 
tant with molecular weight sa 1526; its critical micellar 
concentration is 2000 ppm at 23°C [l9|. The kinematic 
viscosity v of the Triton X-305 solution is approximately 
that of pure water, 0.89 centistokes, and its surface ten- 
sion (7d ranges from that of distilled water (73 mN/m) at 
the lowest <f> down to 50 mN/m at <f> = 1000 PDMS 
has a density p = 0.968 g/cm and the samples used had 
surface tensions er^ = 2.25 to 20.0 mN/m and kinematic 
viscosities v — 10 to 1000 centistokes. 

We visualize the spreading droplets using shadowgra- 
phy [13] through gel substrates prepared in clear plastic 
petri dishes. A parallel beam of light passes through the 
material from below; an image is formed on a ground 
glass above the layer and captured by a digital camera 
and framegrabber. In addition, to visualize the arm to- 
pography we extract a calibrated cross section using light 
scattering in gels containing 1 /im polystyrene spheres 
[2ll ]. Each droplet is ejected from a micropipette held 
about 1 cm above the surface of the gel. We examined 
droplets with volume V = 0.1 to 10 /iL. For reference, 
a V = 5 [iL of fluid with a molecular size of 10 -8 m 
(PDMS), corresponds to a monolayer surface film with 
radius ~ 12 cm. 

For the starbursts, the key physics involves a compe- 
tition between spreading due to surface tension and the 
gel rigidity, ultimately resulting in the fracture of the gel 
substrate. A droplet on a gel surface will spread if it 
has a spreading parameter a g — Od — <J g d > 0, for gel- 
air, droplet-air, and gel-droplet surface tensions respec- 
tively. Since both the gel and the droplet are largely 
composed of water, we assume that a g d is negligible and 
define Act = a g — Ud- The forces that drive the ra- 
dial spreading of the drop deform the gel substrate, as 
shown schematically in the inset to Fig. [2l and we es- 
timate the energies involved. Suppose that there is a 
nascent crack of width d (corresponding roughly to the 
length of a stretched polymer at the point of failure) 
in the radial direction, and depth h into the gel. This 
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FIG. 2: Mean number of arms on individual droplets as a 
function of (a) droplet surface tension a&\ (b) gel elastic 
modulus G; and (c) MX (see Eq. [T}, assuming Xpdms = 
\2\Triton- Inset: Schematic of Marangoni/elastic instabil- 
ity. Open symbols are Triton X-305 solutions, with a d = 72 
mN/m (<), 71 mN/m (o), 70.6 mN/m (*), 70 mN/m (>), 
63 mN/m (□), 60 mN/m (x), 59 mN/m (y), 57 mN/m (o), 
55 mN/m (*), 53 mN/m (A), 50 mN/m (+). Solid symbols 
(here overlapping) are PDMS: a d = 20.1 mN/m and v = 10 
cS (■), a d = 20.8 mN/m and v = 50 cS (♦), a d = 20.9 mN/m 
and v = 100 cS (•), a d = 21.2 mN/m and v = 1000 cS (*)• 
Symbol colors represent gel substrate: 1.3 Pa (blue), 3.7 Pa 
(magenta), 5.6 Pa (cyan), 8.2 Pa (black), 11.7 Pa (yellow), 
16 Pa (green), 28 Pa (red). Droplet volume is 2.5 fj,L. Each 
point is averaged over at least ten drops at the same param- 
eters. Horizontal error bars are due to G and a g ; vertical 
bars are the standard error for the measurements. No data 
at (j> — 1000 ppm is included, due to the unreliable counting 
of numerous arms. 



crack extends in the plane of the gel over some fraction 
of the outer circumference of the drop, aR. There are 
now several relevant energies: the energy associated with 
the surface tension stretching E a ~ AcraRd, the elastic 
energy at fracture, E e ~ GaRdh, and the thermal energy, 
Et = ksT. Here, we assume that the failure modulus 
of the gel is proportional to the elastic modulus, G. The 
time to nucleate a crack, r nuc ~ exp^—E^/Ex), involves 
an energy barrier depending on E ei E a , and Et- As- 
suming that d is comparable in size to an extended gel 
polymer strand, and that h is also microscopically small, 
the first two energies are both many orders of magnitude 
larger than Et- Consequently, we expect (and observe) 
rapid nucleation. During r ntIC , fluid from the drop dif- 
fuses a distance A ~ (DTnuc) 1 ^ 2 into the gel. A then sets 
the depth for the crack: h = A. Therefore, we might ex- 
pect a Marangoni/elastic instability above a critical value 
of 
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We characterize the Marangoni/elastic instability in 
terms of several simple measures: the number of arms 
that form, their mean width and their length as a func- 
tion of time. Fig. [2] indicates that the parameter, M cap- 
tures a key feature of the experiments, namely the num- 
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ber of arms. This number increases with Act/G; here we 
omit the factor A, although we expect that it depends on 
the diffusivity of the droplet in the gel. However, data for 
such diffusivities is not available. The number of arms in- 
creases for decreasing G (connect identical symbols, best 
seen for the solid stars and solid line), or for decreas- 
ing ad (for example, dashed lines connecting open sym- 
bols of the same color). In the former case, onset occurs 
by changing the gel concentration for a constant PDMS 
droplet, and in the later by changing the surfactant con- 
centration of the Triton X-305 droplet for a constant gel. 
There is scatter in the plot due to a number of uncon- 
trolled factors, for instance heterogeneities in the gel and 
the impact velocity of the droplet. Remarkably, in spite 
of the three orders of magnitude difference in viscosity 
for the two types of droplets, the only adjustment re- 
quired to achieve coincidence with the Triton X-305 data 
was to assume Xpdms — 12Axriton- This suggests that 
PDMS diffuses more slowly than Triton X-305, in accor- 
dance with its longer molecular length. An alternative 
framework in which to view the instability is to consider 
that the number of arms in the Marangoni/elastic in- 
stability decreases with increased surface tension of the 
droplet. This is to be expected since the wavelength of 
the instability increases with surface tension in general 
for interfacial phenomena flal ]. 

Several limits of Eq. Q] are instructive. For droplets 
with sufficiently high surface tension (low concentration 
of Triton X-305, low Act), the central droplet does not 
spread significantly and the gel remains intact, result- 
ing in no arm growth. Similarly no arms appear for gels 
with G > 70 Pa; instead, the droplets spread as on a 
solid 0, Ell- For 30 Pa < G < 70 Pa, a different mor- 
phology involving the growth of thin, branching wisps 
was observed for Triton X-305 droplets (but not PDMS), 
as shown in Fig. [JJd and 14| . Conversely, for gels with 
G < 1, the substrate is essentially a viscous fluid and 
the droplet spreads out rapidly, as on a deep liquid layer. 
In Fig. [21 note that the open blue symbols (performed 
on gels with G = 1.3 Pa) have a consistently low num- 
ber of arms, suggesting that the mechanism is nearing 
its limit of applicability: the gel no longer has a yield 
strength. For the remainder of this work, we will focus 
on the starbursts. 

For each of the major arms in a starburst pattern, we 
track the length L(t) as a function of the time from its 
initiation, as shown in Fig. [3l After ~ 10 seconds, the 
evolution of the pattern is substantially accomplished, 
and the scaling discussed below typically breaks down. 
For each arm, we fit L(t) to L = At a . Fig. 2] shows 
that the fit value of the exponent a is consistently close 
to 3/4, independent of substrate gel modulus, droplet 
surface tension, arm width, and volume. In addition, 
the factor of 1000 in viscosity does not appear to affect 
a, which is indistinguishable for the PDMS and Triton 
X-305 data. 
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FIG. 3: Length of the five arms (A,B,C,D,E) of a Triton X- 
305 starburst (inset) at a function of time. Time t is measured 
from the start of each arm's growth, not the drop release time. 
Solid lines are fits with a free exponent; dashed lines are fixed 
a — 3/4. Gel is 3.7 Pa and droplet has V = 5^L and a d = 60 
mN/m. 
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FIG. 4: Fit value of spreading exponent a vs. (a) shear mod- 
ulus G, for constant V = 5fiL, (b) surface tension o^, for 
constant V = 5fiL, (c) arm width, for constant V = 5/iL, and 
(d) droplet volume, V, for constant gel and droplet proper- 
ties. Points are averaged over the arms of a single droplet. 
Legend as given in Fig. [2] 



This time dependence is the same as for oil spreading 
axisymmetrically on water where the value of the 
exponent arises due to the competing effects of viscous 
drag and surface tension forces. In the liquid-on-liquid 
case, the prefactor 



coc(act) 1 / 2 (p M r 1/4 



(2) 



has a weak dependence on the viscosity fi (as well as 
density p, which we did not vary significantly in our ex- 
periments). Intriguingly, recent work on the surfactant- 
induced fracture of a particle raft [l2| , saw a multi-armed 
morphology similar to the starbursts observed here and 
^3/4 S p rea( Ji n g behavior, which the authors associated 
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FIG. 5: Mean spreading prefactor C as a function of mean 
arm width. Each point is averaged over the arms of a single 
drop, with max/min values shown as vertical bars. See legend 
in Fig.[2]for parameters of each point. Size of symbol indicates 
droplet size, with the majority taken at V = 5/iL. 



which the arms advance. The initially smooth gel cracks 
under the tangential stress from the surface tension gra- 
dient imposed by the droplet. Once this process has be- 
gun, the crack opens further due to both the curvature at 
the tip and surface tension forces due to that curvature. 
The tip velocity is set by the viscous spreading of the 
droplet fluid into the opening crack, with resulting dy- 
namics similar to that of the liquid-on- liquid case (Eq.[2]). 
Time-resolved measurements of the initial stages of star- 
burst formation will be necessary to further clarify these 
issues. The novelty of this cracking instability is that it 
is driven by an advancing layer of fluid. 
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Andreotti, and L. Mahadevan for helpful conversations. 



with liquid-on-liquid behavior. 

In order to obtain more precise values of C, we as- 
sume that a — 3/4 is universally applicable. We then 
determine the time-averaged value of the prefactor using 
C = (L(t)/t 3 / 4 ) t ; this parameter sets the scale of the 
spreading rate. Its value is approximately the same for 
both PDMS and Triton X-305 droplets in spite of the dif- 
ference in viscosity, echoing Eq. O Although there does 
not appear to be any significant dependence of C on vis- 
cosity, the average value of C is linearly related to the 
average arm width for a given droplet. We note that 
the arm width at half- length increases and/or decreases 
somewhat during arm growth, again likely due to hetero- 
geneities in the gel. Since the variation in width over the 
course of the growth is small compared to the length of 
the arms, we consider its time-averaged value. As shown 
in Fig. [5l the wider the arms of a droplet, the faster they 
grow (larger C). Remarkably, data for the various gel 
moduli, droplet volumes, surface tensions, and droplet 
material (PDMS vs. Triton X-305) are all consistent with 
this single linear relationship. 

In addition to instabilities due to tangential stresses 
described in Eq. [TJ normal stresses due to gravity may 
be relevant as well. We consider the ratio 



pVg 

(Ac) R 



(3) 



and find a value of about 5, indicating that the gravi- 
tational and surface tension effects are similar in magni- 
tude. Although the number of starburst arms slowly in- 
creased with droplet volume, we observed starbursts for 
the smallest volumes that we could generate (see small 
starburst at left side of Fig. [TJ:). 

A key issue is the nature of the mechanism through 
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